Reissner-Nordstrom Black Holes in Extended Palatini Theories 
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We study static, spherically symmetric solutions with an electric field in an extension of general 
relativity (GR) containing a Ricci-squared term and formulated in the Palatini formalism. We 
find that all the solutions present a central core whose area is proportional to the Planck area 
times the number of charges. Far from the core, curvature invariants quickly tend to those of 
the usual Reissner-Nordstrom solution, though the structure of horizons may be different. In fact, 
besides the structures found in the Reissner-Nordstrom solution of GR, we find black hole solutions 
with just one nondegenerate horizon (Schwarzschild-like) , and nonsingular black holes and naked 
cores. The charge-to-mass ratio of the nonsingular solutions implies that the core matter density is 
independent of the specific amounts of charge and mass and of order the Planck density. We discuss 
the physical implications of these results for astrophysical and microscopic black holes, construct 
the Penrose diagrams of some illustrative cases, and show that the maximal analytical extension of 
the nonsingular solutions implies a bounce of the radial coordinate. 
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I. INTRODUCTION 

According to General Relativity (GR), the fate of any 
sufficiently massive star is to end up its lifetime form- 
ing a black hole, a spacetime region that contains a 
zero- volume singularity of infinite density cloaked by an 
event horizon. Singularity and uniqueness theorems, to- 
gether with the cosmic censorship conjecture tell us 
that black holes are described by solely three parameters: 
mass, charge and angular momentum, a result known as 
the no-hair theorem that yields the Kerr-Newman fam- 
ily For any black hole, the event horizon acts as a 
sort of no- way-out layer that separates the physics out- 
side the horizon, which is the one that we can directly 
explore and find it in excellent agreement with the GR 
predictions Q, from the physics of its interior, where 
much less is known. When quantum phenomena come 
into play, the mere existence of an event horizon induces 
the emission of thermal particles via Hawking radiation 
[H • Quantum effects may also affect the details of gravi- 
tational collapse, as recently studied in Q , and are likely 
to play a very important role when the spacetime cur- 
vature reaches the Planck scale. In fact, it is generally 
accepted that classical black hole singularities should be 
removed by quantum gravitational effects. However, our 
current understanding of quantum gravity is not mature 
enough to provide a clear and detailed description of how 
this could occur. It seems thus justified to explore phe- 
nomenological approaches to this problem hoping to gain 
some new insights that help us to better understand how 



black hole structure could be modified by quantum grav- 
ity effects. 

As the singularity theorems only state some precise 
physical conditions under which the appearance of sin- 
gularities is unavoidable, several approaches have been 
carried out to find conditions that avoid or ameliorate 
the formation of black hole singularities. A well known 
example is Bardeen's black hole @ , in which exotic mat- 
ter sources get rid of the singularity while keeping the 
horizons and the asymptotically flat character. This sin- 
gularity avoidance is realized through the formation of a 
central matter (de Sitter) core, such that the correspond- 
ing spacetime is interpreted as the gravitational field of a 
nonlinear magnetic monopole and can be derived from a 
nonlinear electrodynamics model Q • That approach be- 
came a prototype for most developments on regular black 
holes within classical GR. In this sense, there has been 
much activity aimed at finding alternative matter sources 
for the interior region, such as introducing nonlinearities 
as in the case of nonlinear theories of electrodynamics 
8], implementing a de Sitter core that matches the ex- 
terior field in some "junction" region (see e.g. [1, [To| 
and references therein), or using new ideas inspired by 



noncommutative geometry 11], in such a way that sin 
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gularities are removed. Also other regular magnetically 
charged solutions within GR have been found 

A different approach comes from the idea that exten- 
sions of GR with high-curvature corrections could be able 
to capture some essential features useful or required to 
find an effective description of the quantum gravity dy- 
namics. This point is supported by the study of quan- 
tum fields in curved spacetimes j 1 31 ] and by approaches 
to quantum gravity based on string theory [141 . Il5j . In 
this sense, gravity theories containing higher-order curva- 
ture invariants naturally appear as modifications of GR 
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in these quantum gravity approaches [13H15I ]. Such the- 
ories generically lead to higher-order partial derivative 
equations, which is a manifestation of the fact that new 
high-energy degrees of freedom arc being taken into ac- 
count. The resulting solutions are thus expected to de- 
pend on a larger number of parameters (integration con- 
stants), which should provide more freedom and/or new 
mechanisms to avoid the singularities. However, nonsin- 
gular black holes of this kind are rather scarce (see e.g. 

Ill El) 

and usually require the addition of exotic sources 
of matter again. This suggests that the addition of new 
degrees of freedom in the problem is unable by itself to 
resolve the problem of singularities. Moreover, these the- 
ories suffer from ghosts and other perturbative instabil- 
ities [3]. These problems, however, can be avoided if 
the curvature invariants appear in appropriate combina- 
tions, because then the equations of motion may remain 
second-order like in GR. These are known as Lovelock 
gravities (see, for instance, [2(| for a pedagogical in- 
troduction) . For example, the simplest extension of Ein- 
stein gravity via higher-curvature terms in this context 
corresponds to Gauss-Bonnet gravity and picks up three 
new terms, R^ p<J R' ll/p ' y - ^R^W + R 2 : which supple- 
ment the Einstein- Hilbert Lagrangian. It should be noted 
that in 3 + 1 dimensions these new terms are topological 
invariants not contributing to the equations of motion, 
which means that Lovelock gravities only provide modi- 
fied dynamics in the context of higher dimensions. Exact 
static spherically symmetric solutions to Gauss-Bonnet 
theory in vacuum [2l| and with electrostatic fields [22j 
are known, but they still contain singularities or are ill 
defined. 

An alternative strategy to obtain modified gravita- 
tional dynamics beyond GR is to assume that the metric 
and affine structures of the theory are independent [20| . 
When the connection is not constrained a priori to be 
given by the Christoffcl symbols of the metric, one finds 
that even f(R) extensions yield second-order field equa- 
tions, which contrasts with the usual (metric) formula- 
tion of those theories [23| and the general belief that only 
Lovelock gravities have second-order equations. This ap- 
proach, known as Palatini formalism, has been recently 
used to obtain a covariant action 24) for the effective 
Hamiltonian dynamics of loop quantum cosmology [25[ , 
an approach to quantum cosmology based on the non- 
perturbative quantization techniques of loop quantum 
gravity In the Palatini approach [27]], metric and 

connection are regarded as independent entities and the 
field equations are obtained by independent variation of 
the action with respect to both of them. Though this 
does not affect the dynamics of GR Q, it does have im- 
portant consequences for extensions of it [27| . In general, 
one finds that in Palatini theories the matter plays an ac- 
tive role in the construction of the independent connec- 
tion, which ends up producing modified dynamics. When 
there is no matter, the field equations boil down to those 
of GR with an effective cosmological constant, which de- 
pends on the form of the particular Lagrangian chosen. 



This property has made these theories very attractive for 
cosmological applications. 

The unusual role played by the matter in the con- 
struction of the geometry in Palatini theories together 
with the second-order character of their field equations 
makes them specially interesting to explore the effects 
of new gravitational physics on the structure of black 
holes. In this sense, in Ref.[28[ we considered Pala- 
tini f(R) modifications of GR in interaction with modi- 
fied matter sources, through nonlinear electrodynamics 
(NEDs). In Palatini f(R) theories, the modified dy- 
namics is due to a number of new terms on the right- 
hand side of the equations that depend on the trace T 
of the energy-momentum tensor of the matter. Unlike 
Maxwell's electrodynamics, the stress-energy tensor of 
NEDs possesses, in general, a nonvanishing trace, which 
makes them specially suitable to excite the Palatini mod- 
ified dynamics in electrovacuum scenarios. In the par- 
ticular case of Born-Infeld NED coupled to the gravity 
theory f(R) = R ± R 2 /Rp, where Rp is the Planck cur- 
vature, we found that the degree of divergence of the 
Kretschmann scalar near the singularity can be much 
weaker than in GR. 

In this work we go beyond Ref.[28[ and explore how 
black hole structure is affected by new physics at the 
Planck scale by considering a Palatini theory of the form 
f(R, Q) = R+ aR 2 /R P + bQ/R P , where Q = R^W*" , 
i? M „ is the (symmetric) Ricci tensor, and a and b are 
constants. Terms of this kind have been considered in 
the metric approach in an attempt to find ghost and 
singularity free theories of gravity [29|, a nd also in the 
study of black holes coupled to NEDs [17| . In the Pala- 
tini framework, f(R,Q) theories yield second-order field 
equations that exactly boil down to the usual Einstein-de 
Sitter equations in vacuum (see details in (30j and below) , 
which guarantees the absence of ghosts and other dynam- 
ical instabilities. The presence of a Ricci-squared term in 
the action is very important because it leads to modified 
dynamics even for traceless matter sources, such as radi- 
ation and the usual Maxwell electromagnetic field, which 
contrasts with f(R) theories. As a result, this quadratic 
f(R,Q) model provides the simplest extension beyond 
GR of the usual (nonrotating) Reissner-Nordstrom black 
hole. The quadratic f(R, Q) Palatini model proposed 
here has already been studied in cosmological scenarios, 
where it was found that the big bang singularity is re- 
placed by a cosmic bounce in isotropic and anisotropic 
(Bianchi I) universes filled with standard sources of mat- 
ter and radiation [3lj . 

In Refs. [HI, Hll we reported on several key aspects of 
nonsingular black hole solutions found in this model. In 
this paper we shall go further on this subject extending 
those results and providing a comprehensive account of 
all the derivations and technical details. The paper is 
organized as follows. In section [TT| we recall some basic 
aspects of Palatini f(R,Q) theories and write the asso- 
ciated field equations. The particular actions for matter 
(Maxwell) and gravity sectors of our theory are intro- 
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duced in section Hill where we construct all the relevant 
geometric objects. In section IIVI we comment on the 
choice of ansatz for the line element, and proceed to solve 
the metric field equations. Section [V] is devoted to the 
analysis of the interior region and to the characteriza- 
tion of the event horizons. In particular we construct 
the Penrose diagrams and their maximal analytical ex- 
tensions associated to the different black holes and naked 
cores found. The physical aspects of these solutions, in- 
cluding singular and nonsingular black holes, and micro- 
scopic naked cores, are discussed in section fVTl We con- 
clude in section I VIII with a summary and some future 
perspectives. 

II. ACTION AND FIELD EQUATIONS 

We define Palatini f(R,Q) theories as follows 

S[g,r,tp m ] = i J d 4 xy/^gf{R,Q) + S m [g,if) m ], (1) 

where n 2 = 8nG, S m [g,ip m ] represents the matter ac- 
tion, g a p is the space-time metric, R = g^ v R^ v , Q = 
gV a g v P R^Ra/3, R^ = R f ' Mpi/ , and 

R a '' pnv — d/xF"p — d v r™p + r" A r^g — r^r^ . (2) 

Variation of ([T]) with respect to metric and connection 
leads to the following equations [34j 

/rR^u — -zg^v + 2/q-R mq -R q ' v — k 2 T^ v (3) 
V/3 [V~g <JngT + VqR^)] = , (4) 

were we have used the short-hand notation fx = dxf- 
In the above derivation we have assumed a symmetric 
Ricci tensor, R^ v = R V a, and vanishing torsion. The 
condition on the Ricci tensor, R[^ v ] — 0, forces the con- 
nection components V av to be the gradient of a scalar 
function, T° v = d u <f>. In the usual formulation of GR, 
where the connection is given by the Christoffel sym- 
bols of the metric, one finds that = <9 M In y/^g. In 
our theory ([I]), the condition Ru v ] — is equivalent 
to assuming that can be solved as the Levi-Civita 
connection of an auxiliary metric h^v, which leads to 
= <9 M In \J—h. The explicit relation between and 
<7 M „ follows from the field equations and will be discussed 
later. A reason to set R\p. v \ = is that then the field 
equations of ([T]) in vacuum boil down exactly to those 
of GR (with possibly a cosmological constant, depending 
on the function f(R,Q) chosen). This guarantees that 
there are no new propagating degrees of freedom and, 
therefore, the resulting theory is not affected by ghosts 
or other dynamical instabilities. In regions containing 
sources, the equations of ((T|) differ from those of GR due 
to the presence of new matter/energy-dependent terms 
induced by the mismatch between and g^ v , which 
leads to modified gravitational dynamics without intro- 
ducing new dynamical degrees of freedom or higher-order 



derivatives of the metric. If the condition R^] = is 
relaxed, then T^j must also have a purely vectorial com- 
ponent, r™ = dfj,(f>(x) + Bft, which adds new dynamical 
degrees of freedom to the theory. In that case, the dy- 
namics of (Q} differs from that of GR even in the absence 
of matter/energy sources (see [35| for related results in 
this direction). 

We now focus on working out a solution for (jj). At 
first sight, since R and R^ v are functions of the connec- 
tion and its first derivatives, Q can be seen as a non- 
linear, second-order partial differential equation for the 
unknown connection. However, there exist algebraic rela- 
tions between R, R^ v and the energy-momentum tensor 
of the matter that make the problem easier to handle. To 
proceed, we first define the matrix P, whose components 
are P^ = R m g ah ', which allows us to express ([3]) as 

/flP/ - (V + 2f Q PfP a v = « 2 V • (5) 
In matrix notation, this equation reads 

VqP 2 + IrP ~~i= k 2 T , (6) 

where T is the matrix representation of . Note that R 
and Q are the trace of P and P 2 , respectively. The solu- 
tion of © implies that P can be expressed as a function 
of the components of the energy-momentum tensor, i.e., 
P = P{T). Assuming that for a given f(R,Q) theory 
such a solution exists, Eq.Q can now be seen as an al- 
gebraic equation for the connection in which, besides the 
metric g^ v , there is an explicit dependence on the energy- 
momentum tensor of the matter. To solve it, we look for a 
metric h such that (g]) becomes V/3 [\Z^hW] = 0. This 
guarantees that the independent connection can be ex- 
pressed as the Levi-Civita connection of h. Using matrix 
notation, we have 

V^hh- 1 = S~gg- 1 ^ , (7) 

where we have defined X = (/rI + . Computing 

the determinant of this expression, we find h = gdctE. 
With this result, we have 

h- 1 = i , h= Ndctt) ±- x g . (8) 
V det £ V ' 

This shows that the connection of f(R, Q) theories can 
be explicitly solved in terms of the physical metric g^ u 
and the matter sources. 

With the above results, the metric field equations can 
be rewritten in a more compact and transparent form. 
Expressing © as 

Pt = Li + n 2 f , (9) 
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and using (|8]), we can rewrite P Al Q E Q , l/ as R^ah^V det E. 
This allows to express © as 



-0„. + K J„ 



VdetE V 2 ^ 



(10) 



where T/ = T^g"" '. 



III. ELECTRICALLY CHARGED f(R, Q) BLACK 
HOLES 



A. Matter Lagrangian 

The Schwarzschild black hole is the most general spher- 
ically symmetric, nonrotating vacuum solution of GR 
and also of ([3]). However, that solution assumes that 
all the matter is concentrated on a point of infinite 
density, which is not consistent with the dynamics of 
((3]). In fact, if one considers the collapsing object as 
described by a perfect fluid that behaves as radiation 
during the last stages of the collapse, explicit compu- 
tation of the scalar Q = R^R^ shows that the en- 
ergy density p is bounded from above by p ma x = Pp/32, 
where pp = 3c 5 / AttHG 2 ~ 10 94 g/cm 3 is Planck's density 
[H 11 . Therefore, one should study the complicated 
process of collapse of a spherical nonrotating object to 
determine how the Schwarzschild metric is modified in 
our theory. For this reason we study instead vacuum 
space-times with an electric field, which possess a nonzero 
stress-energy tensor able to excite the Palatini dynamics 
even in static settings. The resulting solutions should 
therefore be seen as Planck-scale modifications of the 
usual Reissner-Nordstrom solution. Let us thus consider 
as the matter source in action (p} the Maxwell lagrangian 



S m [9,i>m] = / d 4 x^g~F a pF a ? (11) 

whose associated stress-energy tensor is written as 



rp V _ 



p a p 1 



F ^Fn a 
r a r fi re 



(12) 



where F^ v = d^A y — d v A^ is the field strength tensor 
of the vector potential A^. For a purely electrostatic 
configuration and taking a spherically symmetry line el- 
ement of the form ds 2 — gudt 2 + g rr dr 2 + r 2 dfl 2 , with 
dQ 2 = d9 2 + sin 2 Qdip 2 , one finds that the only nonvan- 
ishing component is F tr . It is then easy to see that the 
(sourceless) field equations V ' ^F^" = lead to 



ptr 



'I 



1 



r 2 \J-gugrr 



(13) 



where q is an integration constant that represents the 
charge generating the electric field. With this result, the 



product Fu a F a u in (fT2"|) becomes 



p ap v _ f-gttg r r(F tr ) 2 I 0\_ ^1 
M a - \ 6 6 ~ [ 6 6, 



(14) 



where / and represent the 2x2 identity and zero ma- 
trices, respectively. Using this, we find that 



-I 



87rr 4 V / 



(15) 



In order to write the field equations in the form (|10l) . 
we first need to find the explicit form of P for this matter 
source, which will allow us to construct E and compute 
its determinant. To do this, we use (fTSj) and write © as 



x 2 _i 6 
6 \ 2 J 



(16) 



where X 2 ± = (/ + ^ ± ^) /2 and we have defined 

k 2 = n 2 /4ir = 2G. It is easy to see that there are 16 
square roots that satisfy this equation, namely, 



/siX- 











S2A- 











s 3 A 


V o 









(17) 

where Si = ±1. Agreement with GR in the low curvature 
regime (where fn — > 1 and /q — > 0) requires Si = +1. 
For this reason, we simplify the notation and take 

From this it follows that the matrix E is given by 
^ /«;. r^fX-I \ (oA 6 \ 



"~ 2 1 1 v ^ jy V 6 x+ij V a + l) ' y ' 

where cr± — + \/2JqX±^ . From this expression it is 

easy to see that V det E = cr + (T_ and, therefore, the field 
equations (pT0|) become 




(20) 



B. Gravity Lagrangian 

To work out the explicit form of the metric we must 
specify an f(R, Q) model. It is very useful to consider 
the family f(R, Q) = f(R) + l 2 p Q, where l P = ^JhGj<? ~ 
10 _35 m is Planck's length, because tracing ^ with the 
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metric g^ v leads to the well known f(R) relation Rfu — 
2/ = k 2 T, which implies that R = R(T). Since for 
Maxwell theory T = 0, it follows that the f{R) part 
of the f(R,Q) theory does not play a very relevant role 
in the dynamics. We will just assume that the f(R) part 
is close to GR, f(R) — R + CL2R 2 + a 3 i? 3 + . . ., and that 
R(T = 0) = for simplicity (and for consistency with 
the choice = +1 above). This is true, in particular, for 
the model 



f(R,Q)=R + l 2 P (aR 2 + Q) 



whose cosmological dynamics has been carefully studied 
in the literature [3l| and that we set as the model to 
be discussed from now on. For this model we have that 
when is given by (fT5j) then R = 0, /r = 1, and 
f(R, Q) = l 2 P Q. Using this in (TTgl) and taking the trace 
we find 



+ l P Q 



k 2 q 2 l 2 



from which we obtain 



Q 



-4 4 

k q 



+ 1%C 



k 2 q 2 lp 



(22) 



(23) 



which coincides with the expression of GR. From this 



result it follows that X± = 



J_ ( I ± ~± 



1 ± ^-V 1 

becomes 



Noting also that / ± 



- 2 2 - 2 2 

K q 1 K q 



and cr± = 



R^(h) 



~i 2 
K q 

2r 4 



±1 6 ' 
6 i/, 



(24) 



which fully determines the dynamics of our f(R, Q) the- 
ory coupled to Maxwell electrodynamics. These equa- 
tions exactly recover GR in the limit lp — > 0. 



IV. SOLVING THE FIELD EQUATIONS 



of radial coordinate r ^ f. This can be seen from the 
relation g^ w = Yi^ a h av / v detS, which leads to 



fgtt 







g rr 

r 2 





r 2 sin 2 ( 





h rr 
<?+ 

















a- / 



From the line element (f2"5j) it 



hgg/<T- implies that 



,-2 



(27) 

is easy to see that ggg = 

k 2 2 l 2 

■ = t 2 — \ F . It is also 
easy to see that g rr = h rr /a + = (h?? / a + )(df / dr) 2 . Since 
the time coordinate is the same in the two line elements, 
for the gu component we just have gu = h t t/<7+- 

It should be noted that, in general, the line el- 
ements can be written without using r and f as 
coordinates [37[ ■ In that case, we would have 
d~s 2 = h ab (x°,x r )dx a dx b + ^(zV 1 )^ 2 and ds 2 = 
g a b(x° ,x 1 )dx a dx b + r 2 (x° , x 1 )dQ 2 , with the relations 
g a b(x°, x 1 ) — /i a f,(x°, x l )/a + and f 2 = r 2 er_, being the 
latter clearly independent of the choice of (x ,^ 1 ). At 
this point, it is important to note that the radial func- 
tion f 2 vanishes when r 4 = k 2 q 2 l 2 P . This means that the 
2— spheres of the ft,— geometry (the geometry associated 
to the independent connection) can only be put into cor- 
respondence with the 2— spheres of the g— geometry up to 
r/(k\q\lp) 1 / 2 > 1. Since the h— geometry deviates from 
the g— geometry by the effects of a matter-induced defor- 
mation (represented by the matrix E^"), the impossibil- 
ity of mapping any portion of the h— geometry into the 
r/(h\q\lp) 1 / 2 < 1 sector of the g— geometry suggests that 
the matter (the electromagnetic field in this case) cannot 
penetrate in that region of the physical spacetime. The 
description and analysis of the properties of the hyper- 
surface rj {k^lp) 1 / 2 — 1 will be one of the main goals of 
this paper. 

Working with the Schwarzschild (canonical) coor- 
dinates of above and introducing the ansatz h u = 
—A(r)e 2 ^^ and hff = 1/A(r), the components of the 
tensor R^Qi) become 



In order to solve for the metric g^, it is more conve- 
nient to solve first for h^ u using (|24j) and then transform 



back to g M „ using the relation g^ — Y 1 ^°'h au / V det £ 
that follows from ©. To do this, it is convenient to de- 
fine a line element associated to the metric using a 
set of Schwarzschild-likc coordinates as follows 



ds 2 = h tt dt 2 + hffdf 2 + f 2 dn 2 



(25) 



This line element is formally identical to that correspond- 
ing to the physical metric g^, 



ds z = gttdi 1 + grrdr 2 + r 2 dH, 2 



(26) 



Rt 



1 



~A 



Rf T = 



1 



Aff 
~T 

2V>f 

. 1 / 

~A 



A- 

) + 2ipff 



A r 2 

-T + 4>f + - 
A r 



Aj . 

— ) + 2ipff 



-r + fo) +--r 

A J r A 



Re 6 = ^ [I - A(l + r^jf) - fAf] 



(28) 



(29) 
(30) 



but their relation is not trivial due to the nonconformal 
relation between the two metrics and the different choice 



From (|24|) it is easy to see that the combination R t — 
Rf r — 0. Using (|2"5)l and (f2"9")l . this combination implies 
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that t/>f = — > %j) = constant, like in GR and f(R) the- 
ories. As usual, this constant can be eliminated by a 
redefinition of the time coordinate, leaving a single func- 
tion (A(f)) to be determined. From (f3T)|) we find that 
A(f) satisfies the following equation 



A - fAp] 



~? 2 



1 



2r 2 



k 2 q 2 lj 



(31) 



Using the ansatz A(f) = 1 — 2M(f)/f and the relation 

~2 2 K, q If 



k 2 q 2 l 2 

r~ = r~ zzr^i turns into 



Mr 



~ 2 2 

kt<T 
4r 2 



(32) 



1/2 

Taking into account that df / dr = a+/ <r_ , the above 
expression becomes 



M r = 



~ 2 2 

K q a + 



Ar 2 a 



1/2 ' 



(33) 



which reduces the problem to solving a first-order dif- 
ferential equation in the variable r. The integration 
constant of this equation can be identified with the 
Schwarzschild mass Mq = rg/2 of the vacuum problem 
(q = 0). We can thus write 2M(r) = r$ + AM to em- 
phasize that it is the function AM which encodes the 
electrostatic contribution to the mass function. In order 
to obtain AM, it is useful to introduce some definitions 
to work with dimensionless variables. We thus define a 
length (squared) associated to the charge, r 2 = k 2 q 2 , and 
a dimensionless radial variable z = r / 'y/r q lp. With this 
notation, the metric g^ v can be expressed as 



A(z) 

9tt — , Qrr 



a-A(z 



S 2 zcrJ 



1/2 



(34) 

where a± = 1 ± 1/z 4 , we used the relation f = ra 1 / 2 and 
defined AM/rg = SiG(z), together with 




62 = 



qt-P 



rs 



(35) 



The (dimensionless) function G(z) has a purely electro- 
static origin and satisfies the following equation 



dG 
dz 



z 4 + l 



Wz 



(36) 



We stress that the two scales of the problem, namely, the 
integration constants q and Mq, have been replaced by 
the dimensionless ratios Si and S2 given in (|3"5)) . 



A. Finding G(z) 

The integration of (|36"|) to obtain the function G(z) 
can be carried out straightforwardly using power series 
expansions in two regions of interest. One is the far limit 
z ^s> 1 and the other is Z ~ 1. 



1. Far limit z>l 

To study this limit, it is useful to express dG/dz as 
follows 



G z 



1 1 

+ 



1 



Z 2 Z S J (1- l/z 4 )!/2 



(37) 



Using the binomial expansion (1 + x) a = X^feLo (jc) fl! *' 
the above expression becomes 

G* = X)(-1)* ( / (^ 4fe + ^ 4fe ) ■ (38) 

k=0 



The integration of this expression is immediate and gives 



^(1 + 4*)^ 



(l + 4fc) 
(5 + 4fc)z 4 



From the first terms of this expansion, G(z) ~ — 1/z — 
3/10z 5 , one can verify that when r ^ lp the GR limit 
(i.e. the standard Reissner-Nordstrom solution) is recov- 
ered 



9tt 



1 



rs 



ril 



q l P 



'Il- 



l's 



2;2 



ril 



q'P 



2r 2 



(40) 
(41) 



where the first-order corrections come from the a± func- 
tions rather than from the expansion of G(z). As it fol- 
lows from (|40| and (|4Tj) . as long as rs ~> lp, the location 
of the external horizon in these black holes is essentially 
the same as in GR. How this picture changes for micro- 
scopic black holes will be discussed later on. 



2. Near limit z — > 1 

To study this limit, we first consider the change of 
variable z 4 = 1 + x, which leads to 



G x = \ 
4 



-1/2, 



a (i + xy s,i + X- 1 ' 2 {i + x) 



-7/4 



(42) 



Using again the binomial expansion, we find the following 
solution 



00 



(Z 4 - l) fc+1 /2 



4^ (fc + 1/2) 
11 



-3/4 
k 



-7/4 
k 



2Vz~^l-— (z-l) 3 / 2 + 0[z-l] 5 / 2 . (43) 
6 

To have agreement with the solution (|3l?)) found before, 
we need to add an integration constant /3 on the right 
hand side of this expansion, which leads to 



G{z) = P + \ Vz 4 -l [/f (z) + f7 (z) 



(44) 
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where f x (z) = 2 F 1 [±,\, §,1 - 
function, and /? ps -1.74804. 

The mass function can then be written as 



is a hypergeometric 



M(z) 



= l + 8i[P 



1 



/#(*) + /*(*) 



(45) 

The series expansion provided here can be used to per- 
form analytical studies of the geometry near z = 1 . For 
that purpose, the computer algebra Mathematica pack- 
age xAct [36| will be very useful. 



V. INTERNAL GEOMETRY AND HORIZONS 

To study the internal geometry of black holes, it is 
convenient to use coordinates in which the metric is well- 
defined even at the event horizons. For line elements of 
the form we are using, ds 2 = —B(r)dt 2 +C(r)dr 2 +r 2 dil 2 , 
it is particularly useful to replace the (t, r) coordinates by 
the so-called Eddington-Finkelstein coordinates (v, r*) 
38] that turn the line element into 



ds 2 = -B{r)dv 2 + 2dvdr* + r 2 (r*)dQ 2 



(46) 



where v = t + x with (dx/dr) 2 = C(r)/B(r), and r = 
r(r*) is such that (dr* jdr) 2 = B(r)C(r) = l/a_. From 
our definitions in ([34 
i 



r q lp 1 and we can also 



we have B{r) — -g*^, C(r) — 



; r c z, where r c 
r c z* (for notational convenience, from now 

instead of 



£± 

(T — A(z) ' 

define r 

on we use the dimensionless variables z and z 
r and r*). 

The line element f|46[) puts forward that the geometry 
is fully characterized by the functions B{z) and z(z*). 
The relation between z and z* can be found by direct 
integration and is given by z* = %Fi [— \, 5, f , tk\ z. For 
z > 1, we have z* ps z — l/6z 3 + . . ., whereas for z — > 1 
we find z* m zl + yJz — \ + -^(z - l) 3 / 2 + . . ., where 
z\ = 0Fr [|]/r [f] « 0.59907. The relation between 2 
and z* is monotonic and invertible in the region z > 1. 
From (|36[) we also see that S(z) is only defined for z > 1. 
As we pointed out above, this is a dynamical consequence 
of the theory, not a coordinate problem. To learn more 
about this point, we need to study the properties of the 
geometry as the region z 1 is approached. 



A. Region z — ¥ 1 

The general expansion of the metric function B(z) of 
P6"| around z«l leads to 



1/ _<v 

2 I 5 2 



1 



25i 



1 - 



yjz - 1 - . 
(z - 1) + . . . (47) 



This expansion shows that the metric component g vv is 
in general divergent as z — > 1. However, it also points 



out the existence of a particular mass-to-charge ratio for 
which the divergence disappears. In fact, if we take 8\ — 
Sf ee -1//3, we get 



B(z) 



1 - 



si 



5d-> 



1 - 



2^ 
3<5 2 



(2-1)- 



(--1) 2 + 



(48) 



where, according to ([33|) . Ji/^ = r q /(2lp) » 1 for 
macroscopic black holes. This result implies that when 
5\ = 5* the line element and the metric components are 
finite everywhere. It is worth noting that this feature 
is quite similar to what is found in some NEDs in GR, 
for which a particular combination between matter and 
black hole parameters give rise to metrics which are finite 
everywhere 39]. In our case, it is the nonlinearity in the 
gravitational action (JlJ, as opposed to the nonlinearity 
in the matter sector of NEDs, which gives rise to this 
effect. It should be noted that even though in GR with 
some NEDs the metric may be finite everywhere, some of 
the curvature invariants always blow up at the center of 
the solutions, according to a well known theorem hold- 
ing for electrostatic fields with Maxwell weak field limit 
|40| . To see what happens in our theory, we consider 
next the behaviour of the Ricci scalar, Ricci squared, 
and Kretschmann scalar of the metric g^ u in the regions 
z — > 1 and also z ^> 1. 



B. Curvature scalars 

Using the exact solution for the metric found above, 
we can compute some relevant curvature invariants to 
extract useful information about the geometry. The sim- 
plest such objects are the Ricci scalar R(g) = R fllJ (g), 
the Ricci squared scalar Q(g) = g^ v g af3 R fla {g)R l/ ij{g) , 
and the Kretschmann scalar K(g) = R a p ll v(g)R a ° IJ, ' J (g), 
where R a p^ v {g) represents the Riemann tensor of the 
metric g a p. In GR, these objects have the following ex- 
act values 



Rgr — , Qgr — -a > Kgr — 



12r? 



24r s r 2 



14r£ 



(49) 

Though the Ricci scalar in GR is zero, because the stress- 
energy tensor of the electromagnetic field is traceless, the 
other two scalars are nonzero and, in fact, diverge as 
r — » 0, which signals the existence of a strong singularity 
at r = 0. 

Since the solution that we found for the metric in our 
theory is given in terms of infinite series expansions, we 
find it useful to compute the above scalars in the two nat- 
ural regimes in which those solutions were found, namely, 
in the region z 3> 1 and in the limit z — > 1. When 
z ^> 1, we find the following expansions (r c = \Jr q lp 
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and z = r/r c ) 



48r| 

r.10 



R(g) 
Q(g) 

K{g) « K GR 



1 - 



.11 



lAArsrl 



(50) 
(51) 
(52) 



It is clear that these results recover GR when r ^> r c 
(or, equivalently, z » 1). However, in the region z — > 1 
the behaviour of those curvature scalars is completely 
different from the z » 1 expansions. When z — > 1, we 
find 



r 2 c R(g) 



-4 



1 

2S~ 2 



mi 

Si 



-O(z-l) 
1 

(z- l) 3 / 2 



(53) 
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1 - 



865 2 526! 



951 



3<5 2 
6(5 2 — 5(5i 



3^( Z - 1)3/2 



5|(z-l) 3 



O(z-l) 



O 



(54) 



(z-lf 



16 



88(5? 645i 



1 



35 2 

2 (2gi - 35 2 ) 

3(52( Z _ 1)3/2 

1 



O (* - 1) 



(55) 



45 2 2 (2 



1 



1 



l 

(^1)2 



From these expansions we see that the curvature scalars 
diverge at z = 1 except for those configurations whose 
charge-to-mass ratio satisfies the condition 8\ = 6* , since 
in that case all of them become finite. This is an impor- 
tant result whose physical consequences will be explored 
in more detail later. For now, it should be noted that the 
avoidance of the singularity is a local, nonperturbative ef- 
fect that has no impact on the structure of the black hole 
at distances z 1, which quickly tends to that of GR 
regardless of the particular value of 6%, as can be seen 
from the expansions in (|50)) . (jSTj) . and (|52)l . 



C. Properties of the hypersurface z = 1 

Consider the normal to a hypersurface 
S(v,r* ,9,<p) ^constant, namely, N = g^ v d v Sd^. 
US = r, then = ^(1,5,0,0), = (0,1,0,0), 

and N^Nn = {^f B = 1/C(r) = a-B(z). This 
result shows that the vector N is spacelike outside the 



external horizon, null at the horizon, and timelike inside 
the horizon except at z — 1, where it becomes null 
again (regardless of the value of Si ) due to the presence 
of <r_. This implies that the singularity found for 
5\ =/= SI lies on a null hypersurface, which contrasts with 
the Schwarzschild (spacelike) and Reissner-Nordstrom 
(timelike) singularities of GR. In the context of GR, 
null singularities have been found in the interior of 
Reissner-Nordstrom black holes perturbed by neutral 
matter [4l|. Further exploration of the connection 
between these two results shall be done elsewhere. 

To learn more about this null hypersurface, we com- 
pute now its surface gravity using the Killing vector 
I — dt- In the coordinates (pl6"|) . which are regular across 
the external horizon, the components of this vector are 
W = (1,0,0,0) and = (-5,1,0,0). Since the surface 
gravity k is defined by V Q (PZ M ) = —2nl a (evaluated at 



z = 1), it follows that 2k — d r *B(r) = a 
leads to 



1/2 



B z /r Cl which 




8<5,*<5 2 r- 



if s x ^ si 
if Si = si 



(56) 



Note that, strictly speaking, the surface gravity only 
makes sense when evaluated on a horizon. Since in the 
cases Si ^ SI the null surface z = 1 is singular, we believe 
that the divergence of k in those cases is just a manifes- 
tation of the breakdown of the geometric description on 
that surface. On the contrary, the vanishing of k at the 
horizon z = 1 when 5± = SI puts forward the smoothness 
of the geometry at that location. We refrain ourselves 
from interpreting these results in a thermodynamic con- 
text because this aspect of Palatini theories of gravity 
has not been considered in the literature with sufficient 
detail yet. 



D. Horizons 

Horizons are located at the points where the metric 
function B(r) vanishes. From the definitions given above, 
this happens when the curves fi(z) = 1 + 5±G(z) and 

f2{z) = S2za 1 / 2 meet. We have already seen analyti- 
cally that for large black holes, r c /r$ -C 1, the external 
horizon lies almost at the same location as in GR. How- 
ever, since the internal structure of our black holes is 
different from that of GR the very existence of an inner 
horizon is not guaranteed a priori. Moreover, the loca- 
tion of the external horizon for microscopic black holes 
may also significantly differ from the prediction of GR. 
For these reasons, in this section we focus on these points 
to complete our analysis of the internal structure of these 
black holes. 

We note that due to the character of the solutions, 
given as infinite power series, the best way to determine 
whether inner horizons exist or not is by using a graphi- 
cal representation of the functions fi(z) and f2(z). From 
these representations, see Figs. |TJ [31 and [31 one realizes 
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that the first terms of the expansions may provide use- 
ful information to understand the main features of the 
various cases of interest. In this sense, it is worth not- 
ing that both fi(z) and fa(z) are monotonic functions 
whose asymptotes are fi(z) ~ 1 — 5\/z and fa{z) ~ S^z, 
respectively. This implies that if for some zo we have 
fa(zo) < fi( z o) then at some zjj > z$ we will have 
fa(zu) = fi{zii), which implies the existence of an hori- 
zon (as can be verified graphically). Following this rea- 
soning, we expand fi(z) and fa{z) around z = 1 and 
identify fa (z) /fa (z) > 1 as the condition for the existence 
of an inner horizon in that region. Using the expansions 
fa{z) « 1 - ft + 25iV7^T and / 2 (z) « 2<5 2 \/F^T, we 
find the following cases: 

• If ft = ft then fi > 1 . 

This relation translates into the condition r q > 21 p. 
Expressing the charge as q — N q e, where e is the 
electron charge and N q the number of charges, we 
can write r q — y / 2a em N q lp, where a em is the elec- 
tromagnetic fine structure constant. With this no- 
tation, the above condition becomes N q > N q = 

y/2/a em ~ 16.55. This result means that in or- 
der to have an horizon when ft = ft*, the num- 
ber of charges must be at least equal or greater 
than N q . For smaller values of the charge one can 
verify graphically that there are neither inner nor 
outer horizons, which represents naked core solu- 
tions. When the charge is greater than N q then we 
have only an external horizon and no inner horizon. 
All of these solutions are free of curvature singular- 
ities, as results from the analysis of section IV Bl 



Horizon location as JV 9 changes when 5i=t>i* 




Figure 1. The location of the external horizon is given by 
the intersection of the curve /1 (solid green) with (solid 
red curves labeled by N q ). The dashed (orange) straight lines 
that meet at the origin correspond to f% R — 82 z. The dashed 
black curve is fi R = 1 — Si/z. Note that /1 and f% quickly re- 
cover the GR behaviour for z>l. As a result, the location of 
the external horizon almost coincides with the GR prediction 
for N q > 30. 

• if ft < a? then g£ <i + M .ff;_V • 



This condition indicates that once the charge-to- 
mass ratio Si < S* and the number of charges N q 
are specified, one can always find some z > 1 that 
verifies the inequality, which implies the existence 
of an horizon. Therefore, regardless of the value 
of N q , when Si < SI we always have an (external) 
horizon (see Fig[2]). 



Horizon location as N q changes when <Si><Si* 




Figure 2. Same notation and labeling as in Fig[T] When Si < 
Si, we have G(z — 1) > 0, which forces all the curves f2(z) 
to cut fi(z) in a single point. Except for very small values 
of N q , the location of the (external) horizon almost coincides 
with the GR prediction. In this plot 5i = — 5 x 1CT 1 ). 

. If Si > 6t then g < 1 - 2St l\-% /2 ■ 

From this it follows that for some combinations of 
Si and N q there may or may not exist a z > 1 that 
satisfies the inequality. This means that in some 
cases there is no horizon (fa(z) > fi(z) always), 
which implies a naked singularity, while in other 
cases there may be up to two horizons. This is 
verified graphically in Fig. [3l where we can appre- 
ciate solutions without horizon, solutions with two 
horizons, and solutions with one (degenerate) hori- 
zon (extreme black hole), a situation analogous to 
that of the standard Reissner-Nordstrom solution 
of GR. 



E. Penrose diagrams and analytical extensions 

The previous discussion about horizons provides useful 
information for the construction of the Penrose diagrams 
corresponding to the solutions found. For instance, we 
have seen that when Si < 6* , all solutions represent black 
holes with a single horizon. Behind the horizon we find 
a singularity located at r = r c . This configuration is es- 
sentially the same as that found in GR for Schwarzschild 
black holes, except for the fact that the singularity is now 
null instead of spacelike. As a result, the corresponding 
Penrose conformal diagram is that represented in Fig. 2] 
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Horizon location as JV 9 changes when 6\ 




Figure 3. Same notation and labeling as in Fig. Q] When 5i > 
S* , some configurations have two horizons (see the curves with 
N q — 50, 75, 100) while others have no horizons (like N q = 
30). The limiting case N q = 35 (not shown here) has only one 
(degenerate) horizon and represents an extreme black hole. 
The set N q < 35 represents naked singularities. In this plot 
Si = <5*(l + 3x 10" 1 ). 





r r 



Figure 5. Penrose diagram for the case Si > 5* with two 
horizons (r + and r~). The only difference with respect to the 
GR diagram is that the singularity lies on a null surface at 
r — r c . 



Figure 4. Penrose diagram for the case <5i < 6*. Unlike in 
the Schwarzschild black hole, where the r — singularity is 
spacelike, the singularity here is null and appears at r = r c . 



When Si > 5*, we may find the same subcases as in 
GR, namely, solutions with two horizons, with one dou- 
ble (degenerate) horizon, and naked singularities. Like 
in the previous example, the main difference is that the 
singularity is null rather timelike. To illustrate how the 
Penrose diagram is modified, we plot the case with two 
horizons in Fig. [5] 

The case Si — 6* deserves special attention because 
the null hypersurface r = r c is nonsingular. Unlike in 
the other cases with Si ^ 5* , the absence of a singularity 



suggests that the geometry may admit some analytical 
extension beyond that point. This extension is naturally 
obtained from the relation {dr* jdr) 2 = l/cr_ that defines 
the function r 2 (r*) in (IH51) . In our analysis following (|46p . 

1 /2 

we implicitly assumed that dr* jdr = l/er_ , and omit- 
ted the possibility of having a branch with the negative 

1/2 

sign, dr* jdr = — l/c_ . In the singular cases Si ^ S*, 
the omission of the branch with dr* jdr < is justified 
because there is no natural way to extend the geometry 
beyond the singularity at r — r c . However, if there is 
no singularity at r = r c , the divergence of dr* jdr at this 
point simply states that the function r(r*) has reached a 
minimum at the point r* — 0.59907r c (see the discussion 
following Eq. (|4"5|) and Fig. [5]) . For values of r* < r* , the 
branch with dr* jdr < describes a new region in which 
the area of the 2— spheres grows as r* — > — oo. The rela- 
tion between r and r* can thus be written explicitly as 
follows (see Fig. 



2F1 
2r* 



1 1 
1F1 



if r* 
if r* 



<r* 



(57) 
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In terms of r* , we also have dG/dz* = a+/z 2 , which 
gives continuity to the metric across the bounce. Using 
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Figure 6. Representation of the curves z = z(z*), dz/dz* and 
dG/dz* . The minimum of 2(2:*) occurs at z* « 0.599, where 
dz/dz* vanishes and dG/dz* reaches its maximum value. Re- 
call that r* — r c z* . 

T\ and r 2 to parameterize the 2— spheres when r* > r* 
and r* < r*, respectively, the line element (|4"6")l can be 
written as 



ds 2 



-B(n)dv 2 



jdvdri + r 2 dil 2 



-B(r 2 )dv 2 Yj2dvdr 2 



if r* 
if r* 



>r' c 
<r* 



This representation is useful to understand that for fu- 
ture directed (dv > 0) worldlines, drdv < on timelike 
or null worldlines if r* > r*, but drdv > if r* < r*. 
From a physical point of view, this means that a spheri- 
cal shell of matter that collapses and crosses the external 
horizon will unavoidably shrink to a sphere of area Airr 2 
before bouncing off into an outgoing shell of increasing 
area. The outgoing shell, obviously, cannot return to the 
region from which it was sent initially because it crossed 
an event horizon. As can be seen from the Penrose dia- 
gram of this spacetime in Fig. the outgoing shell can 
reach several different final regions (different / + regions 
for light rays). 



Figure 7. Penrose diagram for the nonsingular case Si =51. 
The null surfaces labeled as r c represent the region where the 
area of the 2— spheres reaches its minimum (nonzero) value. 
An incoming null geodesic (light ray) propagating from Im 
towards the horizon r + , will reach r c and bounce off as an 
outgoing null geodesic into another region, eventually reach- 
ing I^ L (if no scattering takes place). The surface r c , there- 
fore, plays the role of (attractive) black hole and (repulsive) 
white hole at the same time. 



admits a nice physical interpretation since it suggests 
that each charge sourcing the electric field has associated 
an elementary quantum of area of magnitude \/2a em Ap. 
From this it follows that the ratio of the total charge 
q by the area of this surface gives a universal constant, 
p q = q/(4nr 2 c ) = (4tt\/2) -1 ^fc 1 /(hG 2 ), which up to a 
factor \pl coincides with the Planck surface charge den- 
sity. This result is independent of the mass of the black 
hole and, therefore, applies both to singular as well as to 
regular solutions. 

If we focus now on the regular solutions only, we find 
that the regularity condition 8\ = SI sets the following 
mass-to-charge relation 



VI. PHYSICAL ASPECTS 

We now discuss several points useful to understand 
some physical aspects of the solutions described in pre- 
vious sections. 




r S = W7Z\ 7- 



{r q l P f /2 



1 mp 
16*^ < (59) 



/ 3 



which can be interpreted in the sense that the matter 
density inside a sphere of radius r c = (rqlp) 1 / 2 becomes 
another universal constant 



A. Quantum nature of the nonsingular solutions 



We first note that using r q = y/2a em i\ q ip 



that the area of the r = 
N q y / 2a em Ap, where Ap 



N q lp we find 
surface is given by A core = 
Airlp is Planck's area. This 



pp 
4(5* 



(60) 



From the definition of 5\ and /o core , comparing any two 
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configurations one can verify that 



(«) 



o {b) 

pcore 

~Ja)~ 

Pcore 



(61) 



This relation is very useful to classify the different black 
hole configurations, because the relation (1 + (35i) that 
appears in (|47|) can be written as (1 — 8i/o~l) = (1 — 
Pcore/ Pcore), where p* core is the density of the regular core. 
This representation allows for a more physical interpre- 
tation and classification of the solutions given in Sec. [V] 
In particular, the nonsingular configuration p core = p* ore 
is the only case in which the core density is a magnitude 
independent of Mq and q and, in fact, turns out to be 
given in terms of fundamental constants only. This fact 
puts forward the special quantum-gravitational nature of 
the nonsingular solutions. In this sense, we believe that 
the regularity condition p core — p* ore , rather than as a 
fine-tuning issue, should be interpreted in the spirit of 
a quantum constraint relating the mass and charge (or 
charges, in general) of the solutions to avoid the singu- 
larity. It should also be noted that adiabatic changes of 
the charge and mass of the black hole [42j do not allow 
transitions between nonsingular configurations, which in- 
dicates that such configurations have different quantum 
numbers (different q and M but the same p* core )- This 
must have important consequences for Hawking radia- 
tion, because if the emission of quanta is to be compatible 
with the regularity of the core, then the spectrum must 
necessarily have a discrete structure. We shall leave the 
exploration of this issue for future works. 



B. Physical and analytical extensions of the 
electrostatic solution 



The above results picture a black hole interior with 
an ultracompact core of radius r c which contains all the 
mass in its interior and all the charge on its surface. This 
view is physically very appealing but should be compat- 
ible with the mathematical solution represented in Fig. 
[7J In fact, Fig. [7] represents an exact mathematical solu- 
tion of a physically incomplete problem, because it does 
not take into account the necessary existence of the mas- 
sive charged particles that generate the electrostatic field. 
The fact that the surface r = r c is null suggests that if 
massive charged particles were explicitly included in the 
problem, then the Killing vector field dt could become 
again timelike in the region hidden by the r = r c hori- 
zon, where the sources should be located, thus allowing 
for the existence of static interior solutions of the type 
suggested by the constraint (|59p . Therefore, for the de- 
scription of the geometry behind the r = r c horizon, one 
should specify the of the sources that carry the mass 
of the core and the charge that generates the external 
electric field, which would allow to define a new auxiliary 
metric h^ v able to parameterize the internal geometry 
of the core (assuming that suitable matching conditions 



can be found at r = r c ). Since the regularity condition 
([55]) supports that the core matter density is bounded, 
we expect the existence of completely regular solutions 
in agreement with the results found for this same theory 
of gravity in cosmological scenarios [3l| . 

It is worth pointing out that if the r = r c null surface 
could be smoothly matched to an interior region filled 
with matter, where dt were timelike, then the analytical 
extension of the exterior spacetime beyond r = r c would 
admit two possible and different branches, which would 
appear on different sheets of a larger conformal diagram. 
One sheet would contain the matter-filled region and an- 
other the analytical extension shown in Fig. where 
the radial coordinate bounces. This could give complete 
physical reality to the spacetime of Fig. [7] in such a way 
that particles approaching r — r c would be scattered into 
the white hole region instead of falling into the matter- 
filled sector. This and related phenomenological issues 
will be explored elsewhere. 



C. Astrophysical black holes 

From the large z expansion in (J30J) and pTj) . we saw 

2 

that the GR solution B GR (r) = I-^ + ttt is an excellent 
approximation for any r 3> lp. This implies that the 
location of the external horizon of these charged black 
holes is essentially the same as in GR. We thus find 



I£(l + ^l_2r»/r|) = (62) 
l + y/l-ASl/(N qy /2 



rs 
2 

i_s_ 

2 



where (|59|) has been used. For a solar mass black 
hole, where the number of protons is around N Pi q ~ 
10 57 , Eq.(|59p implies that the number of charges needed 
to avoid the z = 1 singularity is just N q> Q = 
{2r s 5l/lp) 2 /*/j2c7^ a 2.91 x 10 2 



«2(i 



~ 484 moles), 
which is a very tiny amount on astrophysical terms. In 
fact, this amount of charge is comparatively so small, 
Nq t& /N P}Q ~ 10~ 31 , that it seems reasonable to ex- 
pect that a quantum gravitational violation of electric 
charge conservation could naturally act to avoid black 
hole singularities in stellar collapse processes. In gen- 
eral, N q = N q ,o(M /Mq) 2 / 3 implies that in astrophysical 
scenarios r + rs (see Figs. HJ [2] and [3] ) . 



D. Microscopic black holes 

We have already seen that as N q drops below the crit- 
ical value Nq = ^fT/ a em f=a 16.55, the external horizon 
disappears and the core undresses becoming directly ob- 
servable (see Fig. [1}. The Penrose diagram correspond- 
ing to the maximal analytical extension of these naked 
regular cores, N q < N°, is depicted in Fig. [51 
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Figure 8. Penrose diagram for the nonsingular case 5i = 5* 
without event horizon, N q < Nq. If the asymptotic region I^ L 
is identified with I^ R , then the surface r c should behave as a 
perfectly reflecting boundary. In the general case, however, 
light rays reaching r c would go through into a new universe. 

From it follows that the mass of these objects is 

Nl / \2a em f/ i m P N* /2 m P 
M= 46* ■ (63) 

For the particular case N q — N°, we find M c = 
m P /(V25l) = TO P 7r 3 / 2 /(3r[3/4] 2 ) « 1.23605m P . It 
is worth mentioning that this very particular number 
also appears in the computation of the total energy as- 
sociated to the electric field in Born-Infeld nonlinear 
electrodynamics, which indeed is found to be e(q) — 
7r 3/2//3pp/ 4 ]2^i/2 g i/2 with p the Born-Infcld param- 
eter [39(. Its presence here puts forward the regulariz- 
ing role played by gravitation. Moreover, as the electro- 
magnetic field itself is expected to receive corrections at 
high energies (see [43[ for some string theory results in 
this regard) , it also points out that an improved descrip- 
tion of the electromagnetic field through nonlinear elec- 
trodynamics in these scenarios could provide interesting 
new insights on the interaction between the nonlineari- 
ties of matter and gravitation. All these facts suggest 
that the theory ([TJ) may shed new light on the prob- 
lem of sources in electrodynamics coupled to gravitation 
[32l HH , which already arises in the Reissner- Nordstrom 
and Schwarzschild solutions of GR [hH ]. 

Finally it is worth noting that if N q is seen as an inte- 
ger number, besides the core area, the mass of these non- 
singular naked cores and black holes is quantized, as has 
been recently claimed in [44| on very general grounds. As 
pointed out before, one thus expects a discrete spectrum 
of Hawking radiation, because physically allowed transi- 



tions should occur only between regular configurations. 
This illustrates how Planck-scale physics may affect the 
perturbative predictions of the semiclassical approach. 
New investigations on all these issues are currently un- 
derway and shall be published elsewhere. 

E. Our model as a deformation of GR 

From the analysis of the geometry in the z — > 1 re- 
gion, we have found that the solutions of our theory can 
be classified into three types according to the value of 
Si or, equivalently, according to the value of the core 
density p CO re- When p core > p* core , then the confor- 
mal diagram of Fig. 0] shows a structure very similar 
to that found for the Schwarzschild black holes of GR. 
This admits a nice physical interpretation because when 
the core contains more mass than in the nonsingular 
case, the black hole looks more like an uncharged ob- 
ject, such as a Schwarzschild black hole. On the other 
hand, when the mass-to-charge ratio is smaller than the 
expected value of a regular configuration, p core < p* core , 
then the resulting structure is closer to a typical Reissner- 
Nordstrom black hole of GR (see Fig. [5J. Our black 
holes, therefore, are somehow smoothly interpolating be- 
tween the abrupt structures found in GR. The nonsingu- 
lar solutions, p core = p* corel lie in the middle of these two 
disconnected branches (Schwarzschild - Vs - Reissner- 
Nordstrom), and represent a kind of object not found 
in GR but that emerges when Planck scale effects are 
incorporated in the problem. 

The smooth interpolation found here between the 
Schwarzschild and Reissner-Nordstrdm solutions of GR 
shares some resemblance with the behaviour of the con- 
tracting and expanding branches of the nonsingular cos- 
mological solutions found in this Palatini f(R, Q) model 
[3l| . In GR, one can classify cosmological solutions in 
two groups, one corresponding to an expanding branch 
that starts in a (big bang) singularity, and its time rever- 
sal, which corresponds to a contracting branch that ends 
in that (big crunch) singularity. In the quadratic Pala- 
tini model studied here, the two singular branches of GR 
are smoothly connected through a bouncing solution that 
avoids the singularity. 



VII. SUMMARY AND CONCLUSIONS 

In this work we have studied the structure of electri- 
cally charged black holes in a Palatini extension of GR 
characterized by a Ricci-squared term. Theories of this 
type could be naturally motivated by quantum effects 
in curved spacetimes. These theories provide modified 
dynamics without introducing new dynamical degrees of 
freedom. This implies, in particular, that the resulting 
solutions can be completely classified using the same pa- 
rameters as one finds in GR, namely, the total charge q 
and total mass Mq. The absence of new dynamical de- 
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grees of freedom also guarantees that these theories are 
free from ghosts and other potential dynamical instabil- 
ities. 

We have obtained exact analytical solutions expressed 
as infinite power series expansions. These solutions show 
that the structure of these black holes coincides with that 
of the well known Reissner-Nordstrdm black ho les of GR 
for values of r ^> r c , where r c — lp^j2N q /N^, with N q 

representing the number of charges and N£ ~ 16.55. Im- 
portant modifications arise as we approach the minimum 
of the radial coordinate, the region r — > r c . At this ra- 
dius, the gauge invariant quantity X = — hF^F^" that 
represents the energy density of the electromagnetic field 
reaches its maximum value X max — ppc 2 /2, and other 
quantities of interest such as T^ u also take (finite) Planck 
scale values. This means that the modified dynamics of 
our Palatini f(R, Q) model has the effect of setting up- 
per bounds at the Planck scale on the energy density of 
the matter fields involved, a property already observed 
in cosmological settings. 

Wc have found that r — r c is a singular null hypersur- 
face if 61 ^ 6% (see (f3"S")) for the definition of Si). When 
the mass-to-char ge r atio <5i is set to the particular value 
51 = 3r[3/4] 2 /v^3 ~ o.572, then this null surface be- 
comes nonsingular and the geometry can be analytically 
extended by means of a bounce of the radial coordinate 
(see Fig. [5]). We pointed out that the area of the null 
hypersurface r — r c grows linearly with the number of 
charges, A core = N q y/2a em Ap. As already mentioned, 
this behaviour suggests that each charge sourcing the 
electric field has associated an elementary quantum of 
area of magnitude \J2a em Ap. Using this result, direct 
computation of the surface charge density, p q = q/(4irr 2 ), 
gives p q = 

(47rv / 2)"V c7 /( ?lG2 )> which U P t0 a factor 
\/2 coincides with the Planck surface charge density. For 
the nonsingular solution, the condition Si — 5$ can be 
seen as indicating that the mass density of the core is 
Pcore = Pp/4#l, i.e., it is of order the Planck mass den- 
sity. It must be noted that p q and p* ore are given in 
terms of the fundamental constants H, G, and the speed 
of light c, and are independent of q and Mq (in the singu- 
lar cases, Pcore does depend on q and Mo). In our opinion, 
this is a clear manifestation of the quantum gravitational 
nature of the nonsingular solutions. Rather than as a 
fine-tuning problem, this constraint on the core density 
should be seen as a quantization condition that selects a 
discrete set among all the classically allowed solutions. In 
a sense, this is analogous to Bohr's atomic model, where 
the stability of Hydrogen under electromagnetic emis- 



sion of radiation was postulated assuming the existence 
of certain privileged orbits that had to satisfy specific 
quantization conditions. 

On the other hand, we have found the mass spectrum 
given in (|63|) . which is valid for all (positive) values of 
N q . This mass spectrum has important implications for 
the emission of Hawking quanta. This is so because if 
physically allowed transitions occur between nonsingular 
configurations only, then the resulting spectrum must be 
discrete. Our analysis also puts forward the existence of 
a new kind of nonsingular objects which are not hidden 
by an external horizon. These naked cores exist for values 
of the charge comprised within the interval < N q < Nq. 

An important lesson that follows from our analysis is 
that the boundedness of the energy density does not nec- 
essarily imply that the spacetime is nonsingular. In order 
to find nonsingular solutions, it is necessary that charge 
and mass satisfy a particular relation. In this sense, it 
seems fair to say that the mass spectrum and sizes of the 
nonsingular objects described here need not be in corre- 
spondence with actual physical particles. From the dis- 
cussion of section IIVI on the relation between r 2 and f 2 , 
it follows that including new scales in the problem (such 
as the masses and different gauge charges of the particles 
making up the system) should have an impact on the 
resulting value of r c , defined as the value of r at which 
the energy density reaches its maximum. The analogous 
of the regularity condition 5± = 5* could also set more 
complicated constraints between the total mass and total 
charges of the system, thus providing a richer structure 
to the set of nonsingular solutions. These aspects, to- 
gether with the process of formation and the stability 
under perturbations of the nonsingular solutions (with 
and without external horizon) studied here is currently 
underway. 
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